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Abstract 



For the orthorhombic intermetallic semiconductor AI2RU, the bandstruc- 
ture, valence charge density, zone center optical phonon frequencies, and Born 
effective charge and electronic dielectric tensors are calculated using varia- 
tional density functional perturbation theory with ab initio pseudopotentials 
and a plane wave basis set. Good agreement is obtained with recent measure- 
ments on polycrystalline samples which showed anomalously strong far IR 
absorption by optical phonons, while analysis of the valence charge density 
shows that the static ionic charges of Al and Ru are negligible. Hybridiza- 
tion is proposed as the single origin both of the semiconducting gap and the 
anomalous Born effective charges. Analogous behavior is expected in related 
compounds such as NiSnZr, PbTe, skutterudites, and Al-transition-metal qua- 
sicrystals. 
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AI2R11 is an intermetallic semiconductor!!] quite close in composition to the Al-transition- 
metal quasicrystals. Because of its relatively simple crystal structure, it has been suggested!! 
that investigation of AI2R.U can contribute to the understanding of the pseudogap in the 
electronic density of states observed in quasicrystals. From first-principles calculations of 
the bandstructure and charge density performed using linear-muffin-tin-orbitals (LMTO) in 
the atomic sphere approximation,! it has been shown that the fundamental gap results from 
the hybridization of sp(Al) and d(Ru) states. More recently, hybridization has been identified 
as the origin of the pseudogap in other aluminum-rich transition-metal compounds.0! 

This high degree of hybridization in states near the fundamental gap might be expected 
to have implications for other properties as well. In perovskite-structure oxides, it has 
been shown that hybridization between oxygen p orbitals and transition-metal d orbitals 
results in values for Born effective charges much larger than the nominal ionic values.i!H 
In AI2RU, experimental measurement of the optical conductivity shows a strikingly large 
oscillator strength for phonon absorption.!'! The estimated values of Born effective charges, 
comparable to those of typical ionic insulators, are very much larger than the static ionic 
charges, which previous calculationi has shown to be negligible. In this paper, we report 
the first-principles calculation of the zone-center phonon frequencies, the Born effective 
charges, and the dielectric tensor of AI2R.U. We obtain good agreement for the far-infrared 
optical conductivity with the experimental measurements. Confirming that the static charge 
transfer between Al and Ru is negligible, we conclude that hybridization may be responsible 
for the large effective charges, providing a single mechanism linking the formation of the 
semiconducting gap and the observed strong phonon optical absorption in the far IR. 

The first principles calculations were performed using the ab initio pseudopotential 
method with a plane wave basis set and the conjugate gradients algorithm.^ For Al, we 
used Hamann-Schluter-Chiang pseudopotentials0 with reference configuration 3s 1 3p 1 3d 1 
and cutoff radii of 1.1, 1.2, and 1.3 a.u. For Ru, we constructed scalar relativistic opti- 
mized pseudopotentials@ with reference configuration 5s 1/i75 5p ' 25 4(i 6 , and cutoff radii of 2.2, 
2.4, and 2.05 a.u. The pseudowavefunctions for s, p, and d orbitals were expanded using 3, 



4, and 6 Bessel functions with g c 's of 4, 4, and 6.05 y/By, respectively. The pseudopotentials 
were put into separable form with two projector sill for each angular momentum with the 
/ = component taken as local for both Al and Ru. The energy cutoff was taken as 500 eV, 
resulting in a kinetic energy convergence of 2 mRy The local density approximation (LDA) 
was used with the exchange-correlation functional of Ceperley and Alder as parametrized 
by Vosko, Wilk and Nusair.0 The conjugate gradients minimization for the self-consistent 
charge density and bandstructure was performed using the program castep 2.1. The force 
constants, Born effective charge tensors Z* a g, and the dielectric tensor eoo were computed 
using the density functional perturbation theory methodlil in the variational formulation.^ 
All calculations were performed for the experimentally observed face-centered orthorhombic 
structure with six atoms per unit cell (space group Fddd), lattice constants a = 8.015 A, 
b = 4.715 A, and c = 8.780 A, and free structural parameter x = |.0 The special k-point set 
was a 4x6x4 Monkhorst-Pack gridEl in the Brillouin zone (BZ) of a simple orthorhombic 
lattice with lattice constants |, |, and |, folded in to obtain 48 k points in the full BZ 
of the primitive face-centered orthorhombic lattice [Fig. 1(a)]. The acoustic sum rule for 
phonons and charge neutrality for Z* aj3 were imposed by adding small corrections (< 1%) 
to the largest matrix elements in each case. 

The first-principles bandstructure of Al 2 Ru is shown in Fig. 1(b), with the main features 
in good agreement with previously calculated bandstructures.iH At the bottom of the va- 
lence band are free-electron-like bands derived mainly from Al s states. In the upper half 
of the valence band, flat bands derived from Ru d states are seen to be strongly hybridized 
with the more dispersive Al s and p bands. These calculations show that Al 2 Ru is a semi- 
conductor, with a 0.35 eV indirect band gap from T to the conduction band minimum along 
r — B. The total valence charge density in the z = plane, shown in Fig. 2(a), is quite 
similar to the valence charge density of superimposed free neutral atoms. The difference, 
shown in Fig. 2(b), can be characterized as a small shift of charge from the core regions into 
the interstitial region between the Al and Ru atoms. From this, it seems clear that there is 
no significant charge transfer between Al and Ru, and therefore, negligible ionic character 



for the bonding in this compound. 

From our first-principles computations of the dielectric tensor, we obtain 



18.9 
22.9 
20.7 



(0.1) 



These large values are characteristic of narrow-gap semiconductors, such as PbTe, with 



32.8.0 From computations of the Born effective charge tensors Z* a/3 = dP a /du 



we 



find for Ru, with site symmetry 222, 



-6.28 
-6.96 
-5.40 



(0.2) 



For Al, in the lower symmetry site 2.., we find 



V 
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3.14 
3.48 ±0.70 (0.3) 
±0.47 2.70 

where the plus/minus signs refer to the two different orientations of the nearest neighbor 
shell for symmetry-related Al atoms. In contrast to the negligible static ionic charges, these 
dynamical effective charges correspond in magnitude to a complete transfer of the fourteen 
valence electrons per unit cell to Ru, leaving the Al ions with an effective charge of ±3. The 
phonon absorption is determined by these dynamical charges, and not by the static ionic 
charges, accounting for the large oscillator strength observed in AI2RU. 

To compare these results quantitatively with the measured optical conductivity, we also 
need the phonon frequencies and eigenvectors at q = 0. The point group is D2/1, with eight 
irreducible representations (irreps)0. The dynamical matrix takes a block diagonal form, 
with a 6x6 block corresponding to pure displacements in the x direction (A and B\ irreps) 
and a 12x12 block corresponding to displacements along y and z (B2 and B% irreps). The 
calculated frequencies of the eighteen phonons are listed in Table I with their symmetry 



labels. There are five optically active phonons in AI2R11, with symmetry labels B\ u , B2 U and 

First, we calculate the real part of the optical conductivity <T\(u) = ^e 2 (c^), using an 
oscillator mo deli 

2 2 
<j\(uj) = -Att E t-7 (0.4) 

where 5j is the oscillator strength of the jth mode 

sj= E ( E Y, M * hz laMv)? (0-5) 

a=x,y,z p=x,y,z K 

and k runs over the six atoms in the unit cell with volume Vq, M k is the ion mass, and 
u)j, e(/y), and jj are the jth dynamical matrix eigenfrequency, eigenvector, and broadening, 
respectively. Since we are comparing with results from a polycrystalline sample, we sum over 
all modes neglecting the polarization index, and divide by 3 to average over all directions. 
The experimental measurement^! at 300 K, shown in Fig. 3(a), includes four peaks. We 
associate these with the five IR-active modes expected from symmetry by assuming that 
the intense peak at 265 cm -1 consists of two unresolved modes. The calculated frequencies 
thus represent underestimates of between 12% and 18% of the experimental frequencies, in 
parentheses: 133 cm -1 (151 cm" 1 ), 216 cm -1 and 248 cm -1 (265 cm -1 ), 279 cm -1 (336 
cm -1 ), and 353 cm -1 (405 cm -1 ). Values of the broadening 71, 72 and 73 for plotting 
the theoretical <Ji(co) are chosen to reproduce the intensities of the lowest two experimental 
peaks, with 74 and 75 chosen to reproduce the qualitative shape of U\ (w) at higher frequencies 
(Fig. 3(b)). While the lowest calculated peak is in excellent agreement with experiment, the 
splitting of the second and third calculated peaks is seen to be slightly too large. In addition, 
their combined oscillator strength, obtained from the calculated eigenvectors, is rather too 
small relative to the upper two first-principles peaks. However, the calculated value for the 
integrated oscillator strength 8 / ai(u)du is 4.74xl0 28 s -2 , in excellent agreement with the 
experimental value of 4.94xl0 28 s -2 . Since this quantity is mainly determined by Z* a g and 
is relatively insensitive to the phonon frequencies and eigenvectors, we conclude that the 
large first-principles values for the effective charges agree quite well with experiment. 



Since the real part of the optical conductivity is independent of the dielectric tensor 
eoo, an experimental value for comparison with our calculations must be obtained from the 
imaginary part of the optical conductivity 02(0;) = —f-ex(u). We fit the measured 02(0;) 
below 100 cm -1 to c\uj + C2UJ 2 . Setting c\ = — 4^00 — ^Sj^f; we obtain an experimental 
value of the polycrystalline average of 17, which is 18% smaller than the average of the 
calculated values \{e xx + t yy + e zz ) = 20.8. This overestimate is typical of the results for this 
quantity calculated using LDA.0 

Our first-principles calculations show that large Born effective charges are an intrinsic 
feature of AI2R.U, consistent with experimental observation of strong phonon absorption. 
However, these large charges are not the result of a static charge transfer typical of ionic 
compounds. Instead, given that the key role of hybridization in opening the semiconduct- 
ing gap in AI2RU has already been established,! we propose that the hybridized nature of 
states near the gap is also responsible for the anomalous effective charges in AI2R.U, in direct 
analogy to the perovskite oxidesBB'i This could be more precisely formulated with a tight- 
binding analysis of the first-principles bandstructure@, though the necessary parametriza- 
tion is complicated by the relatively low symmetry of Al 2 Ru. On a more positive note, 
we expect that similar behavior can be observed in other semiconductors with narrow gaps 
opened by hybridization. Indeed, experimental indications of large phonon absorption are 
available for systems including NiSnZr,@ PbTe,il and skutterudites,i! and theoretical sup- 
port for the non-ionic nature of the effective charges in PbTe and the other rocksalt IV- VI 
compounds has been obtained with the empirical pseudopotential method.il Building on 
first-principles investigations and tight-binding parametrizations of the bandstructures of 
these systemsJlBi! a general understanding of the mechanisms producing anomalous ef- 
fective charges could be achieved. With a local formulation transferable to nonperiodic 
systems, this approach could eventually promote the interpretation of optical conductivity 
measurementsi to elucidate the nature of vibrational modes in quasicrystals. 
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TABLES 

TABLE I. The calculated phonon frequencies at q = 0, in cm , grouped according to the 



symmetry labels given in Ref. 21, The oscillator strengths Sj for the optically active modes are 
given in parentheses. 
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216 (2.2) 


133 (0.84) 
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353 (0.67) 
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FIGURES 

FIG. 1. (a) Irreducible wedge of the face-centered orthorhombic Brillouin zone of AI2RU, shown 
embedded in the Brillouin zone of the simple orthorhombic lattice with lattice constants §, §, § used 
in the construction of the k-point sampling set. (b) Bandstructure of AI2RU along high-symmetry 
lines of the face-centered orthorhombic Brillouin zone shown in (a). The indirect gap of 0.35 eV is 
from the top of the valence band at T to the bottom of the conduction band along T — B. 

FIG. 2. (a) Total valence charge density of AI2RU in the z = plane. The central Ru atom 
is surrounded by a hexagon of Al atoms. Contour intervals are 1.3xl0 -4 electrons per unit cell. 
Counting from Al towards Ru, the contour labels increase from 1 to 32, then decrease again to zero 
at the Ru nucleus, (b) Difference of the valence charge density (a) from that of superimposed free 
neutral atoms (Al s 2 p and Ru cfs 1 ). Contour intervals are 1/3 of those in (a), with dashed lines 
indicating negative contour values. Counting outward from Al, the first few contour labels are -2, 
-2, -1, 0, 1 and 2. 

FIG. 3. The real part of the optical conductivity of AI2RU, in (Q cm) -1 : (a) experimental 
measurement at 300 K (Ref. |9|); (b) theoretical prediction. In (b), a background of 50 (O cm) -1 
has been added and values of jj= 25, 9.0, 4.5, 65 and 50 cm -1 (j = 1, ...,5) chosen as described in 
the text. 
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Ogut and Rabe, Figure 1 



Ogut and Rabe, Figure 2(a) 




Ogut and Rabe, Figure 2(b) 
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Ogut and Rabe, Figure 3 



